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ABSTRACT 

The  dc  and  microwave  responses  of  the  BaxSri-x(X,Y)yTii.yC)3  family  of  ferroelectric 
compounds  with  various  substitutional  additives  X3+,  Ys+  are  analyzed  by  combining  the 
random- field  technique  with  the  mean-field  (Landau-Devonshire)  theory  of  ferroelectricity, 
along  with  a self-consistent  computation  of  the  dielectric  constant  of  the  host  material  in  the 
presence  of  the  impurity  fields.  The  fields  in  the  material  are  assumed  to  arise  from  charge 
compensation  at  the  Ti4+  sites,  leading  to  permanent  dipoles  made  up  of  the  resulting  positive 
and  negative  ions  separated  by  a few  lattice  constants.  It  is  shown  that  whereas  completely 
random  placement  of  positive  and  negative  ions  generates  a Holtsmark  distribution  of  electric 
field,  with  infinite  second  moment  and  hence  extremely  large  fluctuations  in  field  strength, 
the  association  of  ionized  impurities  into  permanent  dipoles  leads  to  much  lower  fluctuations 
in  field  and  a distribution  with  finite  second  moment,  which  makes  a self-consistent  dielectric 
constant  meaningful. 


INTRODUCTION 


The  usefulness  of  ferroelectric  materials  in  the  design  of  electronic  devices  arises 
from  the  ease  with  which  their  electrical  characteristics,  e.g.,  conductivity,  dielectric 
constant,  optical  birefringence,  etc.,  can  be  controllably  modified  so  as  to  produce  a 
desired  functionality.  A standard  way  of  controllably  modifying  a material  is  to  introduce 
foreign  elements  into  a host  material  at  concentrations  that  are  too  low  to  change  its 
intrinsic  chemical  nature,  but  which  modify  its  behavior  at  the  macroscopic  level.  In 
ferroelectrics,  these  modifications  are  usually  associated  with  the  presence  of  built-in 
electric  fields,  which  are  generated  by  thermal  ionization  of  the  added  impurity  atoms. 
This  type  of  processing  is  well  known  in  the  semiconductor  industry,  where  materials  are 
routinely  “doped”  to  make  them  conductive  and  otherwise  electrically  active. 


However,  because  semiconductors  already  contain  large  amounts  of  free,  i.e., 
mobile,  charge,  it  is  normally  assumed  that  there  are  no  macroscopic  fields  inside  them. 
Because  doping  such  materials  with  impurities  creates  ionized  donor  and  acceptor  centers 
in  addition  to  the  free  charge,  it  would  seem  that  long-range  fields  should  appear. 
However,  the  free  charge  can  move  to  screen  out  any  such  fields  [1]  over  distances  that 
are  larger  than  the  Debye  length  in  the  bulk: 


lD  ~ 


cm  at  room  temperature. 


(1) 
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where  n is  the  carrier  density,  e the  dielectric  constant.  For  n~  1012  cm'3 , which  is  a 
lower  limit  on  how  low  the  free  charge  can  be  in  most  semiconductors,  this  gives  /D  = 1.2 
pm,  i.e.,  very  small.  Over  these  distances  the  fields  can  only  create  small  potential 
fluctuations,  which  modulate  the  bottom  of  the  material  conduction  band  to  produce  so- 
called  "bandtails",  as  shown  long  ago  by  Evan  Kane  [2]  and  Halpcrin  and  Lax  [3],  In 
contrast,  insulating  materials  can  have  vanishingly  small  amounts  of  free  charge  in  their 
interiors,  and  hence  can  sustain  long-range  fields.  Although  these  fields  will  vary 
randomly  in  magnitude  and  direction,  they  can  superimpose  and  become  quite  large.  In 
this  paper  we  develop  a statistical  method  to  calculate  the  behavior  of  fields  associated 
with  charged  impurities  in  ferroelectrics,  and  use  it  to  discuss  the  effect  these  fields  have 
on  the  host  material's  dielectric  response. 

FIELD  PROBABILITY  DISTRIBUTIONS 

The  appropriate  description  of  random  electric  fields  in  an  insulator  is  statistical, 
since  the  added  charge  centers  are  located  randomly  throughout  the  material  bulk.  This 
means  that  at  any  given  point  in  the  material  wc  can  assume  that  there  is  a probability 
distribution  for  the  field,  i.e.,  that  the  field  is  a random  variable.  The  usefulness  of  this 
picture  depends  on  what  observable  we  plan  to  calculate.  If  the  observables  of  interest 
involve  volume  averages,  these  can  be  evaluated  conveniently  by  invoking  the  ergodic 
hypothesis,  i.e.,  we  can  replace  the  volume  average  by  an  ensemble  average.  To  proceed 
beyond  this  point  the  probability  distribution  of  the  quantity  being  averaged  is  required. 

We  begin  by  calculating  the  probability  distribution  for  electric  fields  generated 
by  one  species  of  added  impurity.  If  the  impurity  centers  arc  uniformly  distributed,  then 
the  probability  pj  (jc)  that  there  is  an  impurity  at  position  .v  is  simply  l/L,  where  V is  the 
host  volume.  To  find  the  field  distribution  we  use  a method  pioneered  by  Markov  [4], 
discussed  first  in  the  astronomy  literature  by  Holtsmark  [5]  and  then  in  the  astrophysics 
literature  by  Chandrasekhar  [6],  and  finally  in  solid  state  by  P.  W.  Anderson  [7]  to 
discuss  electron  localization  in  disordered  media  and  by  Stoncham  [8]  to  discuss  optical 
line  shapes.  First,  we  note  that  if  the  impurity  positions  jq,^,^,...  are  random  vector 

variables  with  probability  distributions  p\  (Jq ),  p2  (x2 ),  py (.<3 ) and 

/ = f(x  | a:j,  *2,  %,...)  is  some  vector  function  of  these  variables  and  the  position  x of  a 

point  in  the  solid,  then  / is  a random  variable  with  a probability  density  given  by  the 
expression 

P(f(x)  = /o)  = (8(/0-/(4m  ^2.-^3.-))  = J II  pn(xn)dx„S(fo  - /(.v|.vi , .v2 , -V3 ,...))  (2) 

n=l 

In  our  specific  case,  / will  be  the  field  due  to  the  random  charges. 

Using  the  exponential  representation  of  the  8- function 

6(L)=  (3) 

-i  (2tt)3 
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we  can  write  the  following  expression  for  the  density  P associated  with  the  overall 
probability  that  the  electric  field  Ej  due  to  the  charged  impurities  has  a value  E at  the 
point  x : 

°°  j3 

P(E,(x)  = E)=V~2N \Y\dxn  | exp[/<7  {E-E, (i|xj ,x2,x^, ...)}] q-  (4) 

" _oo  (2lt) 


Uncorrelated  Positive  and  Negative  Monopoles 


Assume  now  that  we  put  an  equal  number  N of  positive  and  negative  charges 
(“donors”  and  “acceptors”  in  semiconductor  language)  into  the  material,  so  that  it  is 
electrically  neutral  (charge-compensated).  Let  these  charges  be  at  locations 
Xq\ , Xq2’^D3  ’••■xA\  >xA2’xA3  with  probability  densities  p jjj  (xDj ),  Paj  (*Aj ) that 


are  uniform  and  uncorrelated  as  asserted  above,  i.e.,  p Dj  (xDj  ) = pAj  (xAj ) = — . 
Then  the  function  Ej  (x\x\,X2,Xt,,...)  has  the  form 


N 


Ej(x\xdi,xd2 *A1'*A2>-)  = V s 

j = 1 x~xDj 


XAj 


N _ 

= £ E(x-xDj)-E(x-xAj) 

j = 1 

x 

where  E(x ) = e — — is  the  electric  field  of  a point  charge  and  e is  the  relative 

permittivity  of  the  host  material.  Then  Eq.  (1)  can  be  rewritten  as  follows: 

P(E{  (x)  = E)  = J expfitf  • E]  d q (v  -1 J exp[  -14  ■ E(x)\k)  (l-1  Jexp[/?  • ^(ylldy)^ 

— oo  (2k) 


(5) 


(6) 


If  the  densities  of  donors  and  acceptors  are  no-n, i = n,  we  can  write  V = n/N.  Using  this 
notation,  we  have  that 


(k  1Jexp[-i^-£(5)]tfcJV  =(l-U  1 J"  {l  — exp[— ■ ^(Je)]}^)^ 
This  lets  us  write  [9] 

~ ,3 

P(Ej(x)  = E)=  fexp[;^-  £]exp[-2nReF(<?)] %r 

(27t)3 

Now,  the  integral 

ReF(cj)  = j{l  - cos[<?  • E(x)^lx 


(J) 

(8) 

(9) 
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reduces  the  function  P to  Holtsmark’s  distribution  (see  Ref.  [5])  with  a total  density  of 
carriers  2 n.  Some  manipulation  gives 

ReF(f/)  = ^(2K)3/2c3/2|<r/|3/2  (10) 

which  implies  a distribution  of  the  form 

°°  j3 

P(E,(x)  = E)  = fexp[p£-c|?|3/2J— (11) 

(2tt)3 

where  o = ^-(27t)3|,2«e3^2  = 8.400«e3^2  is  a scale  factor  with  dimensions  [field]1 2. 
This  reduces  to  the  1 D integral  of  Holtsmark 

P(E/(x)  = E)  = — — Jexpf— ct-y3 ^ 2 /|z?|'  [vsin  xdx  (12) 

2tc2 |c|  o 

which  cannot  be  reduced  further. 


If  the  host  material  has  a lattice  constant  a,  we  can  write  n = ca  3,  where  c is  the 

( „ VV2 

impurity  concentration  (dimensionless).  Then  a = 8.400r 


\ea2  J 


introduces  a scaling 


electric  field  Esc  =— into  the  problem,  in  terms  of  which  a = 8.400c£3/2.  Taking  a 


£a 


= 4 A and  £=  1000,  which  are  typical  values  for  ferroelectrics,  we  find  that  this  field  is 
about  1 10  V/pm,  i.e.,  close  to  breakdown  of  the  material. 


Fully  dipole-correlated  charge  distributions 


Suppose  now  that  every  donor  is  spatially  paired  with  an  acceptor,  i.e.,  is 
physically  separated  from  it  by  only  a few  lattice  constants.  The  resulting  pairs  of  charges 
can  then  be  treated  as  unit  sources  of  random  electric  field,  i.  e.,  dipoles.  At  low  carrier 
concentrations  (so  that  the  dipoles  are  far  apart)  we  can  treat  the  fields  of  these  dipoles  as 

- OX 

we  did  the  fields  above.  Now,  however,  the  field  E(x)  = e is  replaced  by  the  field 

elf  I 


of  a dipole,  i.e.,  a field  of  the  form 


E(x)=p 


3x(xn)-x2n 


(13) 


where  n gives  the  direction  of  the  dipole  and  pits  strength.  For  a cubic  crystal  the 
dipole  direction  can  vary  in  6 cubic  directions  for  nearest  neighbors  and  8 for  next  nearest 
neighbors.  Let  us  ignore  the  discreteness  of  the  possible  orientations  and  suppose  that 
pis  a constant  vector  that  can  point  in  any  direction,  i.e.,  the  Heisenberg  problem.  Then 
we  can  replace  the  sum  by  an  angular  average: 
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(14) 


V 1 Jexp[-zV/  ■ E(x)]f£c  — > V 1 -^-jdnjexp[-iq  ■ E(x)]dx 


--V1  f <&—[ 

J 4ttJ 


exp 


-i-^rq  ■ J(x)  dh=  f < S > dx 

Del3  r J 


where 

J ='ih{h-q)-q . 


(15) 


Considerable  manipulation  eventually  gives  the  result 


°°  j 3 

P(EI(x)  = E)  = | expfi#  • E - «y)^|] 


where  y = 


it  2pr 


1 + 


sinh  1 -\/3 
2x/3 


(2it) 


> 5.89p . Unlike  the  Holtsmark  distribution,  this 


(16) 


expression  can  be  evaluated  in  closed  form,  leading  eventually  to  the  result 


P(E,(x)  = E)  = 


y n 


*2(£2+YV)2 


(17) 


which  appears  in  the  work  of  Ma  [10].  Again,  a scaling  field  enters  into  the  problem:  if 
we  take  the  dipole  moment  to  equal  ea,  where  a is  the  lattice  constant,  this  field  is  once 
£ 

more  Esc  = — — , where  now  we  write  yn  = 5.89 cEsc . 
za2 
P 


Figure  1.  Holtsmark  and  dipole  distributions.  The  dipole  curve  is  scaled  by  a factor  of 
1/100. 
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Figure  1 shows  a numerical  plot  of  this  function,  together  with  the  Holtsmark 
distribution,  for  the  same  lattice  constant  and  scaling  field.  Note  that  both  these  functions 
are  centered  around  zero,  since  they  are  vector  distributions  and  there  is  no  preferred 
direction  in  the  crystal.  For  convenience  purposes  the  functions  plotted  are  actually 
9 9 9 9 E 

2tc  a Pff(x)  and  2 n a PD(x) , where  x = , since  the  Holtsmark  function 

a2'3 

2 2 

2k  a Pfj  (jf)  is  a universal  function  in  this  coordinate.  In  contrast,  the  dipole  function 

2 2 2r 

2k  o Pn  ( x ) takes  the  form  — , where  F is  a function  only  of  the 

concentration:  T = 1.43  cl/3.  The  dipole  plot  shown  is  for  c = .05,  on  a scale  reduced  by  a 
factor  of  100.  Note  that  the  Holtsmark  distribution  is  much  broader  than  the  dipole 
distribution,  with  a longer  tail.  This  difference  is  a manifestation  of  the  long  range  of  the 
Coulomb  interaction,  and  in  fact  causes  the  Holtsmark  distribution  to  have  no  second 
moment  - i.e.,  the  fluctuations  are  so  large  that  their  mean  square  deviation  from  the 
average  field  is  infinite!  In  contrast,  the  dipole  distribution  has  finite  fluctuations, 
although  it,  too,  is  somewhat  badly  behaved,  having  no  moments  higher  than  the  second. 

THERMODYNAMIC  VARIABLES 

Let  us  characterize  the  ferroelectric  by  using  the  Landau-Devonshire  theory  [11], 
In  this  picture,  the  properties  of  the  ferroelectric  near  the  Curie  temperature  arc 
determined  by  its  free  energy,  which  is  of  the  form 

F =-ED  + ~a(T)D2  +~bD4  + ~cD6  (18) 

2 4 6 

where  D is  the  displacement,  E denotes  an  externa!  electric  field,  and  a,  b , and  c are 
material  parameters.  We  will  assume  here  that  (1)  the  material  undergoes  a second-order 
phase  transition  at  Tc , which  requires  that  c = 0,  b > 0,  and  (2)  that  only  a depends  on 
temperature,  specifically  that  a(T)  = OXT  - Tc ) , where  Tc  is  the  Curie  temperature.  Then 
the  relation  between  the  displacement  D and  the  electric  field  £ is  found  by  minimizing 
the  free  energy  with  respect  to  the  vector  D : 

^-  = -Ei+[a{T)  + bD2\)i=  0 (19) 

Normally  we  would  just  solve  this  equation  for  D , but  here  the  field  £ is  a random 
variable  so  that  a different  approach  is  needed.  If  we  assume  that  the  crystal  is  isotropic, 
then  the  vectors  D and  £ are  always  parallel  and  D = Dh  where  h is  the  direction  of 
the  field  £ . Then 

D = ]—-E=>D  = ^LeU !—  £ (20) 

a + bD 2 a + bD2'  1 a + bD2 

where  E is  the  magnitude  of  E . This  implies  that 
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■ =(a  + bD2J 


(21) 


In  the  presence  of  an  applied  field  Ea , the  average  displacement  is  given  by 


E f 

< Df  >=  J Df  (E)P(E -Ea)d3E  = -j-SDi(E) 


d3E 


(22) 


Ei  +\E-E^ 


where  E . = 7M  = 5.89| 


' e N 


c is  the  “width”  of  the  field  fluctuation  distribution.  Note 


Vea2  J 

that  the  distribution  is  now  centered  around  the  deterministic  external  field  E„ . 


Dielectric  constant 

Let  us  introduce  the  variable  x = D 2 . Then  E 2 = x(a  + bx)2  and  D = — - — E . 

a+bx 


If  the  z - axis  points  along  the  applied  field  and  we  write  d2E  = 2nE2dE<£, , where 
E,  = cos0  depends  on  the  angle  between  the  integration  field  and  the  applied  field,  the 
integral  becomes 


E f °°  £ 0 * 

< D >=<  h„-D>=  -4-  f — ^—2itEzdE  f 


¥% 


Jt2  Qa  + bx 


~'{e2j  +E2+E2-2EEa% 


2Ef  7 E3dE  ] 

~!—rr  I 


¥% 


(23) 


7t  Qa  + bx_xU2 


Ej  +E2+E2 -2EE 


4 


For  small  external  fields  we  can  expand  the  angular  integrand  in  a power  series.  Then  the 
angular  integrations  are  trivial,  and  we  end  up  with 


2 Et 


<D>= 


E3dE 

8 E 

a + bx 

X£/t£2)3 

E„  + 


8£[-5£2  + 3E2 


5^£2  +£2j 


E3  + ■ 


(24) 


= eEa+VEl+- 

where  e and  (3  are  the  dielectric  constant  and  tunability  of  the  material.  Near  the  Curie 

point  we  can  set  a = 0,  which  simplifies  the  relation  between  E and  x to  E - bx2^2  . This 
in  turn  simplifies  the  integrals,  and  we  find  that 
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(25) 


e = ii£-2/3fc-l/3 
27  / 


P = _Jl£-8%-l/3 
729  / 


Self-consistency 


In  the  spirit  of  mean-field  theory  we  now  specify  that  the  dielectric  constant  used 
to  determine  the  field  fluctuation  width  E ^ be  the  same  as  the  one  we  calculate  from  the 


macroscopic  theory,  i.e.,  E ^ = 5.89| 


f \ 
e 


\c  . At  the  Curie  temperature,  however,  we 


\Ea2  J 

know  from  the  previous  section  that  e - if , from  which  we  obtain 

( 2 V/3 

2l\Ejb\ 


fE.\W 


7 


1 1 36  ec 

— — — . Solving  this  for  E ^ gives  E ^ = 1466| 


r v3 

ec 


\a2  J 


b . Then 


*{TC)=  5.89 


*27 


1 


= .004 


7 


\a2  J 


1 _ 9 

— , i.e.,  the  peak  in  £ is  °c  c . Since  the 


fluctuation  width  is  affected  by  the  external  field  Ea , calculating  the  tunability  requires 
some  care.  If  e - e0  (e  /•)+p(e/')e  2 , we  apply  self-consistency  in  the  following  form: 


e(e f)Ef=  5.89^ -^-1  ^ C0  =>  z(Ef  )=  ^ 


\a2  J 


(26) 


2 2 

Writing  the  width  as  £y  = £yg  + \E£  and  equating  powers  of  E^  gives  a zero-order 
equation 


eo(£/o)-T^ 


70 


and  a first-order  equation 


X = - 


[eMIeM, 


P^/oV/o 

co+£a'[Efo)EfO 


If  e = e0[Ef)+V{Ef 


then 


7 0 


f°'  243  / 


— 8 

At  low  concentrations  we  find  that  the  tunability  is  « c , i.e.,  a very  strong  function  of  c. 


(27) 


(28) 


(29) 
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Numerical  Results 


In  order  to  deal  with  finite  temperatures  we  must  allow  a to  be  nonzero,  which 
complicates  the  analysis.  Making  the  integration  variable  change  E 2 = x(a  + bx )2  and 

D = — - — E leads  to  the  following  expression  for  the  average  displacement: 
a + bx 


2E  °° 

< D >=<  h ■ D >=  — — J ! " (a  + bxf  (a  + 3bx)dx 

-rr 


](e2s  + E2  + x(a  + bx)2  - 2Eax{ l2[a  + bx\ 


(30a) 


The  angular  integral  can  be  carried  out  easily,  but  it  leads  to  an  awkward  expression. 
Again  it  is  clear  that  < D > vanishes  as  Ea  — > 0.  It  is  possible  to  modify  the  theory 
slightly  to  obtain  the  behavior  of  the  material  below  the  Curie  temperature  as  well,  from 
the  integral 


<D>= — ix3l2(-\a\  + bxf(-\a\  + 3bx)dx 

n HA 

} M (30b) 

-1(£/  +£a  +x(-\a\  + bx)2  -2EaxXI2[-\a\  + bx^ 


Figure  2.  Spontaneous  displacement  versus  temperature  for  BaSrTiC>3  with  +3  and  +5 
ions  added.  Impurity  concentrations:  (a)  0.0001,  ( b ) 0.001,  (c)  0.01,  and  (d)  0.1. 
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Numerical  evaluation  of  < D > leads  to  the  curves  shown  in  Fig.  2.  We  find  the 
low-field  dielectric  constant  by  linearizing  in  Eu,  which  yields: 


£ _ 8£/  J.v3  /2(a  + far)3  (a  + 3 bx)dx 
37C  0 [t'j  + \(a  + bx) 2 J3 
Ef  y x3/2(-\a\  + b.xf(-\a\  + 3hx)dx 
3U  \a\/h  [e*  +.v(-|tf|  + fev)2j3 

8 


T>TC  (31a) 


T <TC  (31b) 


Figure  3.  Dielectric  constants  vs.  temperature  for  dipole-doped  BaSrTiOj. 
Concentrations:  (a)  0.005  , (b)  0.006  , (c)  .007  , (<7)  .008. 
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Self-consistency  was  imposed  numerically  on  these  functions  to  obtain  the  curves  shown 
in  Fig.  3.  The  results  are  within  an  order  of  magnitude  of  what  we  typically  observed  in 
experiments.  However,  these  theoretical  curves  do  not  predict  the  experimentally 
observed  shift  in  the  Curie  temperature  with  concentration  [12], 

CONCLUSIONS 

In  light  of  the  crudeness  of  this  theory,  it  is  unlikely  that  it  will  be  very  predictive 
when  applied  to  real  (BatSri-j)(X,Y)Ti03  systems,  for  a number  of  reasons:  the  actual 
phase  transitions  in  these  materials  are  first  order,  they  are  ceramics,  etc.  Including  these 
phenomena  will  require  at  a minimum  the  use  of  a version  of  the  Landau-Devonshire 
theory  with  c * 0,  making  the  analysis  more  complex. 


Logio  [C(Y,Ta)  / C(B-site)] 

Figure  4.  Peak  dielectric  constant  vs.  concentration.  A least-squares  fit  gives 
Log10  e(Tc ) = 2.1  - ,92Log(-^-)  . 

cB-site 
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As  an  example  of  the  disagreement  between  theory  and  experiment.  Fig.  4 shows 
the  concentration  dependence  of  the  dielectric  constant  peak  for  members  of  the  material 
set  (BaxSri.x)(Y,Ta)yTi|.y03.  In  contrast  to  the  slope  of -2  predicted  by  theory,  a slope 
close  to  -1  is  observed  [13],  It  is  worth  noting  that  all  these  materials  had  different  grain 
sizes,  and  that  the  Curie  temperature  was  different  for  every  data  point  shown.  This 
suggests  that  the  correct  dependence  will  only  be  derived  when  the  theory  can  predict 
concentration-induced  changes  in  7c  as  well.  We  hope  to  explore  this  issue  in  a future 
publication. 
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